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We have found an exact expression for the decoherence rate of a Josephson charge qubit coupled 
to fluctuating background charges. At low temperatures T the decoherence rate F is linear in T 
while at high temperatures it saturates in agreement with a known classical solution which, however, 
reached at surprisingly high T. In contrast to the classical picture, impurity states spread in a wide 
interval of energies T) may essentially contribute to F. 
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Solid state nanoscale devices provide one of the most 
promising routes to implementing a scalable set of con- 
trollable two-state quantum systems, qubits, based either 
on spin degrees of freedom of electrons in quantum dots 
Q or on discrete charge quantum states in Josephson 
junctions j^. An unavoidable coupling of each qubit to 
the environment leads to decoherence. The loss of coher- 
ence before a sufficient amount of quantum operations 
was performed would be the major impediment in using 
solid-state qubits in quantum computations. 

Recent experiments on Josephson-junction (JJ) single 
qubits 0, 0] have demonstrated the possibility of per- 
forming hundreds of "quantum operations" (i.e. coherent 
oscillations between the qubit states) before environmen- 
tal decoherence sets in. It is believed, however, that tens 
of thousands of such operations are required for quan- 
tum computation to become a reality [5j so that much 
longer decoherence times should be achieved experimen- 
tally. It necessitates a better theoretical understanding 
of realistic mechanisms of decoherence. 

The most conventional way to describe environmental 
decoherence is based on the spin-boson models (see for 
reviews 0) where the qubit interacts with environment 
represented as a set of harmonic oscillators with a given 
frequency spectrum. An alternative approach is based on 
identifying the dominant mechanism of decoherence in a 
real experimental situation and formulating the adequate 
model. The results obtained in such a model are not 
necessarily reducible to those in the spin-boson models. 
This turns out to be the case for decoherence in charge 
JJ qubits addressed in this Letter. 

It is widely believed that in charge qubits the main 
contribution to decoherence comes from an inevitable 
coupling to charge degrees of freedom, which is also re- 
sponsible for 1// noise observed in such devices 0, 
As charge impurities are spatially frozen as experimen- 
tal temperatures T ~ 30 50mK 0, 0, I3i the most 
probable source of dynamical electromagnetic fields is 
impurities recharging, e.g. due to the hybridization of 
their electronic states with conduction electrons (in the 
metallic electrodes, etc.) The appropriate model simi- 
lar to the conventional model of the spectral diffusion in 
glasses '■^\ is known as the spin-fluctuator model. It has 



already been used 0, 0, 0, 0| for a classical ( "high- 
temperature" ) description of decoherence and noise in 
charge JJ qubits. However, in such a description con- 
tributions from impurities with energy levels outside a 
narrow (of order T) strip around the Fermi level of con- 
duction electrons are exponentially suppressed. Had this 
been the case, decoherence from fluctuating charge impu- 
rities would hardly be seen in experiments as it is rather 
unlikely to find such a fluctuator coupled to the qubit 
in the energy strip so narrow compared to the typical 
(atomic) scale over which impurity levels are distributed. 

In this Letter we develop a complete quantum mechan- 
ical description of decoherence due to charge fluctuators 
and obtain asymptotically exact expression for the de- 
coherence rate T{T) given by Eq. H12|l . In the classi- 
cal "high-T" regime it goes over to the previously ob- 
tained result [II EH while in the low-T regime r(r) 
decreases linearly with T, Eq. I|15|l . having a nontrivial 
non-monotonic dependence on the coupling strength g. 
In particular, the exponential suppression of the contri- 
butions of energetically remote impurities turns out to be 
an artefact of the classical description. The hybridization 
with the conduction electrons responsible for the dynam- 
ical recharging leads also to a quantum broadening of the 
impurity levels which results in the contribution of (in- 
evitably present) energetically remote impurities being 
suppressed only as power-law rather than exponentially, 
thus making it detectable and eventually dominant. 

We consider a model where a charge qubit is coupled 
to impurities with charge states fluctuating due to hy- 
bridization with the conduction band: 
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Here luq is the energy split and Ej is the Josephson cou- 
pling of the two qubit states separated from the higher 
states by the Coulomb blockade energy Ec ^ T; both oja 
and Ej are tuneable which allows one to perform quan- 
tum operations on the qubit. Each localized impurity 
state is characterized by its coupling to the qubit Vi (due 
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to the dipole interaction), its energy Si (counted from 
the Fermi energy ep oi the conduction electrons), and its 
switching rate 7.; — 2T:i>\ti\^ {v is the density of states at 
£f in the conduction band and 1^^^ = Vol -ItikP)- AH the 
parameters e^, 7^ and Vi, vary independently in a wide 
interval of energies. 

Solving the Heisenberg equation of motion for the 
full density matrix with separable initial conditions, 
/5b+q(0) = p(0) ® /5b, where p{t) is the reduced density 



matrix of the qubit and — e 



is the equilib- 



rium density matrix of the bath, one writes p(t) in the 
standard formal representation (see, e.g., 11] ): 
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For the charge qubit under consideration, the opera- 
tional condition is Ej <C |3- We are interested in 
decoherence only so that we restrict considerations to 
the 'pure dephasing' regime 0|, Ej = 0. In the per- 
turbative region, small Ej would not lead to noticeable 
corrections to decoherence. Although Ej could be tuned 
to a large value if an operation is performed on the qubit, 
this should happen only for a short period of time which 
gives only negligible corrections to decoherence. 

For Ej — the coupling ct^V" the diagonal elements 
of p do not evolve while the time evolution of the off- 
diagonal elements of p can be represented as 
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where (. . .)b is the average with the Gibbs density ma- 
trix of the bath, pB- Expressions of this sort can be 
exactly calculated in certain problems, e.g. the orthog- 
onality catastrophe 0| and full counting statistics 
in ID, in techniques which are always problem-specific. 
In present considerations, we employ the linked-cluster 
expansion within the Keldysh formalism (similar to that 
used in [l^ for bosonic environment) to find the deco- 
herence rate F defined by 



F(r) = - lim t-^\n\D{t)\ 
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In order to calculate F, we represent D{t) in Eq. Q as the 
following functional integral over the Grassmann fields 
defined on the standard Keldysh contour [l^: 



rVdVdVcDc 
D{t) = I exp 

where the action densities are given by 
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The field Vj{t') = ±Vj for < t' < t, with '-h' sign on the 
upper and ' — ' sign on the lower branch of the Keldysh 
contour, and vanishes for i' < or i' > t. The normal- 
ization Z is defined as the same functional integral but 
with Vj = 0. The integration over the conduction elec- 
tron fields Ck, Ck reduces the action Q to the impurity 
term Sij with the mass operator 
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where gk{t',t") is the conduction electron Green func- 
tion obeying the equation {idt' — £k)5k(^',^") — 5{t' ,t") 
with the delta function defined on the Keldysh contour. 
Now the integration over the fields di and di reduces 
D[t), Eq. ||2Jl, to the appropriate matrix determinant thus 
yielding the following formal result for F, Eq. Q): 
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where G obeys {idf — e° — G = / , and Tr implies an 
integration over times along the Keldysh contour and a 
summation over indices labelling the fluctuators. The de- 
pendence on the running time t above is via Vj (t) defined 
after Eq. (0). The long-t limit in Eq. can be found 
by expanding Trln into power series in iiG. The n*^ or- 
der term of the expansion is a multiple integral over n 
time variables, each running along that part of where 
V ^ 0. As usual 01, we represent each integrand via the 
Keldysh matrices 



f _ I ^ Lt 
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thus reducing each integration over time to that from 
to t (with V — > vTx, where is the Pauli matrix in the 
Keldysh space). Due to the time translation invariance, 
each G depends only on the difference of its time argu- 
ments. Then the n^^ order integrand depends on n — 1 
differences in times, while the integration over the last 
time variable produces the overall factor of t. The re- 
maining integrals in time can be extended to the entire 
axis as all the Green functions exponentially decay in 
time with the time constants 7"^ (or T^^ for the real 
part of the Keldysh component G^). Then the integral 
has a convolution structure in time and, upon a Fourier 
transform, it finally reduces to the integral oft [i?)C/(ijj)] 
over Lo. Since the coefficients of the expansion were not 
affected by the Fourier transform, the re-summing the 
series restores the logarithm and Eq. O reduces to 
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where the trace of In 



1 + ItJ{uj) 



in the space of ma- 



trices © has been explicitly taken and tr refers only to 
the fluctuator matrix indices. 
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It follows from Eq. ® that Tq = r(r = 0) vanishes 
as expected. Indeed, at T = one uses G^{^) — 
[G^{uj) - G^{u;)] sgnw to find 
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since the first order of the expansion is imaginary, while 
the higher orders vanish upon the integration as all the 
poles are in the upper (lower) uj half-plane. 

A further simplification is possible for a typical situ- 
ation when distances between the fluctuators are larger 
than the Fermi wavelength k^^ so that the hybridiza- 
tion is local, iki = tiV~^^'^e'-^''\ Then Eq. © reduces 
to = t^t*g{ri-'Vj,t' — t"). Since g(r) (the Fourier 
transform of the conduction electron Green function g\^) 
oscillates at /cp ^ , off-diagonal matrix elements of E vanish 
upon the integration (neglecting small interference cor- 
rections) making G in Eq. Q diagonal in the fluctuator 
indices, G = GjSij (such a diagonal form was assumed in 
[iQ.] by choosing independent conduction bands for each 
impurity as in 18].) Here 
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23mSr^^- = 27ri/|tjf and ej = e° + SfteE^j is 



the hybridization. The Keldysh component of Q is given 
by H G^(w) = [G^iuj) - Gf{uj)][l - 2n^{uj)l where 



where 7j 

the fluctuator energy (counted from e-p) renormalized by 
; hvbridi: 

ri-piuj) is the Fermi factor with uj counted from e^. 

Then F in Eq. reduces to a sum of the individual 
fluctuator contributions, F(T) = J2j ^ji^)- Substitut- 
ing Gj of Eq. ifTTl) into Eq. ©, subtracting the identi- 
cally zero expression (llUf) for Fo (to improve the integral 
convergency) and taking the real part of the resulting ex- 
pression, we obtain the following contribution of a single 
fluctuator at energy ej = e to the decoherence rate: 
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Here we suppressed the index j and introduced dimen- 
sionless coupling g = v/j of the qubit and the fluctuator 
with the dimensionless density of states Xe{uj) = Trviyeiuj) 
broadened around the energy e by the hybridization: 
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The temperature dependence of the decoherence rate 
(|12|l is presented in Fig. 1 for the strong {g = 100) and 
weak (f/ = 0.2) coupling to a single fluctuator with Ve{^) 
centered at £ = 2e/7 = 3 (or e = 0, insert). In the high-T 
limit, F saturates at the following e-independent value: 



F(oo) = J [1-0(1-5)71^-2 
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FIG. 1: Dependence of the decoherence rate on a tempera- 
ture, Eq. 11211 . for strong and week coupling: the main picture 
is for a fluctuator with v^(lo) centered at e = 3, the insert - 
for £ = 0; here T, F, e are measured in units of 7/2 . 



Although it coincides with the results 0, of the 
classical "high-T" description, it follows from the ex- 
act expression (|12|) only at T ^ max{|e±|,7}, where 
e± = e ± (1/2)\A^2"^— 7^. Thus, for a fluctuator with 
either relatively large e or large v (and thus |e±|), the 
classical regime (|14|l is never reached. Then F is de- 
scribed at any temperature by the "low-T" asymptotics 



F(r) = - arctan^' ^'^ 



TT \ £2 - g2 + 1 



(15) 



which follows for any e from Eq. H12|l either for T <C 
7 or for T <IC min{|e±|}. This means that the result 
of the classical description 0, ^ is applicable only 
for a fluctuator with ^^(u;) centered near the Fermi level 
'provided that T 3> max(w,7). 

Note that a crossover between the asymptotics H15|l 
and p4|l is relatively sharp when A = ma.x{g,^ 3> 1: 
F(T) changes exponentially fast. 



F(r) ^Fooexp(-v47/2r), 



(16) 



in a logarithmically narrow interval. A/ In A < T/^ < A. 

Although a linear in T behaviour similar to that in 
Eq. H15(l would also follow from the spin-boson models 
with the ohmic spectral function, only a full quantum 
treatment of a microscopic model, like that in Eq. , can 
result in a nontrivial T-dependence depicted in Fig. 1. 

One specific and surprising feature of the model (0) is a 
non-monotonic dependence of the decoherence rate F on 
the coupling strength g: at low T a contribution of weakly 
coupled fluctuators can be orders of magnitude higher 
than that of strongly coupled ones, as seen in Fig. 1. Such 
a non-monotonic dependence is depicted in Fig. 2 for the 
fluctuator centered at e = 3. At any flnite temperature 
F as a function of g has a maximum with a cusp at go = 
(1 -he^)^/2_ Only at very high T {f ^ 100 in Fig. 2) 
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tuator must be not far from metallic electrodes to be hy- 
bridized with conduction electrons. Secondly, the peaks 
at e± should be within a few 7's around the Fermi energy. 
This brings further geometrical restrictions essentially re- 
ducing the number of potentially relevant defects so that 
only relatively few fluctuators are likely to contribute to 
decoherence in a typical experimental setup 3]. 
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FIG. 2: Non-monotonic dependence of the decoherence rate 
on the coupling strength at different temperatures for e = 3 . 



the cusp is smeared out and F practically saturates at 
F(oo) — 7/2 as in the classical limit, Eq. (|14|l . 

This surprising suppression is due to an effective spit 
in energy of fluctuators strongly coupled to the qubit. If 
the qubit were in one of its eigenstates, the fluctuator 
energy in the absence of the hybridization would be split 
as e±w/2, Eq. QJ. Allowing for the hybridization, such a 
split to the energies £± = e ± {l/2)^Jv'^ — 7^ occurs only 
for g > 1, i.e. when the coupling v exceeds the width 7 
of the hybridized fluctuator Eq. Ijl^l) . Then the decoher- 
ence of the qubit in the mixed state, Eq. ||2J), is effectively 
contributed from two peaks e± (indeed, the integrand in 
Eq. (|12l) has the two peaks broadened by the hybridiza- 
tion at Lu — e±, besides the exponentially narrow peak at 
Lu ^ 0). Thus, F increases with g at e = until v reaches 
7; a further increase in g pushes the peaks away from the 
Fermi energy, suppressing the hybridization and thus the 
switching rate i.e. effectively freezing the charge states. 
For the arbitrary e, the maximum in F is reached when 
one of the peaks at e± coincides with the Fermi energy. 

The results described by Eqs. 1)12(1 - (|15|l and the sub- 
sequent discussion refer to the case when there are only a 
few fluctuators, so that they can be considered separately 
and their decoherence rates could be simply added. If the 
fluctuators are dense, one needs to average over the rel- 
evant parameters, which is their energies, coupling con- 
stants and switching rates. The averaging, e.g., over the 
energy positions £j, spread within an interval i?, would 
lead to F(T) oc T as long as T <^ E. The reason is that 
the number of effective high-T fluctuators, whose contri- 
bution is described by Eq. (|14|l . would be proportional 
to T while the contribution of each of low-T fluctuators 
would be Hnear in T, Eq. lfTB|l . 

However, as the decoherence rates due to individual 
fluctuators are hugely spread, as illustrated in Fig. 1, the 
effective fluctuators are hardly dense as requirements for 
the effectiveness are rather restrictive. Firstly, the fluc- 
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